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Abstract 



We show in this note how many electron irreducible representations 

of the Lorentz group L can be expressed in terms of Slater determinants. 

v^ , In this way the full configuration wave function of Quantum Chemistry 

is derived without any appeal to the Schrodinger equation and Born- 

f~^ I Oppenheimer approximation. In particular the number nof exchange- 

^SJ . correlations between the orbiting electrons and those in the neutrons may 

^^ ' be calculated and therefore the ground state electron exchange-correlation 

0^ , energy E readily obtained by dividing the atomic energies by n. Com- 

O^ ■ parison with measured atomic energies yields values of E very close to 

C/3 ' -1 hartrees for the 3 neutral atoms Be, B and C. This tallies within a 

O , few percent with the best available E estimates. Although calculations 

C/3 ' were performed for 3 atoms the method is perfectly general and can be 

^~>| applied to any atom. However because E is estimated for metals only, 

wG . gases are not considered in this paper. The intimate relationship with 

t^' density function theory is discussed in the conclusion. 
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^ ! 1 Introduction 

Many years ago the many electron problem was modelled by irreducible repre- 
sentations of the Pauli ring found by constructing the minimal left ideals ^. 
Then when a field was introduced by fibration the well known spin angular 
momentum matrices appeared. The minimal left ideals are the bra vectors of a 
density operator "0 Q , so the method is in fact a density function formulation as 
well as a representation of the two-dimensional unitary group C/2 and therefore 
of a subgroup of the proper orthochronous homogeneous Lorentz group Lp. In 
this paper the many electron case will be extended to the nucleus by considering 
the neutron to be a proton plus an electron so a neutral atom is envisaged as 



an ensemble of A protons and A electrons distributed between N neutrons and 
A-N orbiting electrons. Under these conditions the complex angular momentum 
matrix is |^. 

^ ^ En ®ip +1N ® Ep (1) 

Where 7p, ^n are (P+1)— and (A^+1)— dimensional Lie operators of SO (3) and 
Ep, Ejsf are (P+1) and {N + 1) unit matrices. However the nett electron spins 
must also be matched with the spin of the nuclear ground state even though the 
spins may be orthogonal and therefore not additive. 

According to this scenario the exchange energy of interaction responsible for 
atomic energies will be between the P outer electrons and the P electrons within 
the nucleus and thus the number of interactions will be of the order of P^ unless 
they occur simultaneously. It will be seen that this is indeed supported by 
experimental evidence. Furhermore a preferred spin orbital i/'e will be found for 
the collective motions of the orbiting electrons. 

In order to actually calculate the number of exchange-correlations associated 
with ipe it is necessary to integrate, by exponentiation, the representation of 
the Lie algebra U2 characterising the electronic spin state matching the nuclear 
ground state. When this is done, by means of a matrix exponential theorem, 
the Slater determinants incorporating exchange correlations emerge naturally. 

To illustrate the method, which is quite general, we will begin with ^Be and then 
in Section 3 go on to ^^B and ^^C which belong to another space because A is 
even. Specifically the bosons (even A) are characterised by EucHdean structure 
while the fermions belong to complex manifolds. 



2 BERYLLIUM-9 

According to de Wet [J the many electron configuration decomposes into the 
sum J2 CaPa of projection operators Pa characterising the spin state [A — A, A] 
and coefficients C\ which are polynomials in 7 of equation (1) and may be 
determined by the recurrence formula: 

{n + 1) Cn+i = 7C„ + {A-{n- 1)) C„_i, n>l,A>{n-l) (2) 

where 

7 = Ci =£:„®7P+7Ar®Pp, Co = l (1) 

is related to the usual spin matrices by 7 = 2is. Thus we find 

71 = 1, C2 = i(7' + ^) (2a) 

n = 2, 6C3 = 7 (7' + (3v4 - 2)) = 73 + 257 



= £^5 ® (7p + 257p) + 372, ® 7p + 377V ® 7p + (7w + 257w) «) £;4 



(2b) 



when A=9. Here C3 defines the spin state | corresponding to 6 electrons with 
positive spin and 3 with negative spin, thus matching the spin of the ground 
state of ^Be nucleus where P=4, N=5. In this way the electronic and nucleon 
spins are linked even though their spins may be orthogonal and therefore not 
additive. 

If we use the diagonal representations 

i7P = {4; 2; 0; -2; -4}, ijN = {5; 3; 1; -1; -3; -5} 

then C3 can easily be obtained from (2b). Half the values appear in Table I; 
the remainder are repeated with a change of sign. 
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If however we consider electrons spinning about the X2- axis then 

{jml I 17 I jm) = [{i -m){j + m + 1)] ^ Smf,m+i 
where i7p is characterised by j = 2, ijN by j = |and 
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X 

-X^ 



Which belongs to a complex space even though the eigenvalues of X are real. 
After symmetrization, by interchanging rows and columns, X is found to have 
the eigenvalues of Table I which label its 15 rows, and we will assume from 
symmetry considerations that X77 corresponds to the eigenvalues 1 1| which is the 
same as the lowest non-zero collective eigenvalue due to 4 spinning protons and 
5 spinning neutrons obtained by the method outlined in [4] . This corresponds to 
identifying the 3 lowest non-zero eigenvalues of ^with electrons in the nucleus 
and the high eigenvalue 10. 5i with the collective motion of the 4 orbital electrons. 
Then when (5) is exponentiated the Slater determinants and principal minors 
governing the number of exchange correlations of the 4 electrons with those in 
the nucleus will emerge without any further assumptions. 

The matrix exponential theorem derived by de Wet [3, 4] says that 
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(5) 



o/^e 



E 
= 0,1,, 



i>^kFk{i\k 



■ cos Xk9 



E 
fc = 1, 



Fk{iXk) 



sin A 1.0 



(6) 



where A^ is the k-th eigenvalue oi fj, — X, and 

Fkifi) = i^(M)/(M' + Ai), F,{ti)Fk{fi) = (6a) 

automatically selects an occupied orbital representing an excited state [k] be- 
cause 

F(/i) = ^(a.2 + l){^l^ + Xl) (m^ + XI) = 0. (6b) 

To see it's physical significance we write 

Fkifi) = m'"~^ + Am'""' + /32m'""' + + /3n-iM (7) 

where fi^ -^ "M"^; M^ ^ +/^^) etc. because (5) belongs to a complex space and 

/5i = Af + Al + + Xk\ + A2 , (Afe = 0) (7a) 

/32 = A?A2 + A^Ai + + A?A2 + + ALiA^ , (A^ = 0) (7b) 



/3„_i = A^Ai Xk" A2 = det Xfe, (Afe = 1) (7„) 

Equation (7n) is a Slater determinant of the possible orbits, excluding the state 
[k], and the remaining equations are principal minors representing clusters of 
correlated electrons |^|. The expansion coefficient in the full configuration wave 
function (6) are the matrix powers [m'^]jj. which depend on the state [k] and the 
spins and rotations associated with k. 

To prove (6), differentiate at 6* = to get 
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''=° fc = l,2,. 





is idempotent so ^ Kk{^J.) is a decompositon of unity and (6) follows. The con- 
nection between idemponency and statistics will be discussed in the conclusion. 

Returning to (6) the sin term is characterised by odd powers of // and belongs 
to the matrix (5). The cosine term wih belong to a Euclidean space defined by 
matrices along the diagonal. 

In the case of the problem under consideration the independent eigenvalues are 

\2 _ r 1. 1. 49. 4411 
■^ - \4'-^'T'^/ 

and the relevant matrix powers are 



(8) 
(8a) 



X-n = |, [X^]^^ = 269.390625, [X^]^^ = 29888.34375, [X^]^^ = 3297590.147461 

so (6) yields the electronic wave function 

1/^77 = ^(-22sinf + 8sin6l-3sin|6'+ 15sin^6l) (9) 



which satisfies the criterion -^ 



Comparing (9) with the full configuration interaction wave function proposed 
by Szabo and Ostlund ([5], Ch. 2) the coefficients of the sine terms are integral 
exchange correlations up to a normalization constant (which at least ensures 
that one eigenvalue of the set assumes the value of unity as required by the 
theorem). 

In particular the spin orbital sin(^6') associated with the 4 orbiting electrons 
would have 15 exchange correlations which would be approximately the case if all 
exchanges between these electrons and 4 in the nucleus occured independently 
implying 16 exchanges. Johnson et al.||g| list -14.6674 hartrees for the measured 
exchange-correlation energy of Be which yields an exchange-correlation energy 
per electron of -14. 6674/15=-0. 9778 hartrees. In the next section (Table II) this 
value will be compared with the best available estimate given by the formula of 
Gell-Mann and Brueckner ||7||. 



3 Boron-10, Carbon-12 



The same method may be used with the difference that because we are now 
dealing with bosons that occupy Euclidean space, equation (6) assumes the 
simpler form |^. 

n 

fc = 0,l,.. 

Fk{^i) = ^i{^i -i){ji-\2) {^i- \n)/{p - Afc), F,{^I)Fk{^i) = o 

and 

XkFk{fi)/fiFk{Xk) (10a) 

is idempotent, so that 
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In the case of boron we shall employ the coefficient 

24C4 = 7^ + 5272 + 240 



= ^5 ® (7I, + 52-fl) + A-/N ® (7p + 137P) + 67^ «> 7p + M"f% + ^^1n) «> 7p + 

(7^ + 527^) 0^5 + 240 (11) 

corresponding to a spin state 1, with the degenerate eigenvalue 

A = {210; 42; 10; 2; -14} (12) 

which in terms of the theorem have to be normaUsed so that 1 eigenvalue is 
unity and all eigenvalues are positive. To do this we add At = 14 and divide by 
A/ = 16 to find the positive canonical form 

/i = (C4 + 14)/16 = {14; 3.5; 1.5; 1; 0} (12a) 

The constant At will translate the spectrum because if AX — \X, then 

{A - \t)X = {X - Xt)X 

And after exponentiation will multiply the wave function by a factor e*^'^. 
On the other hand A/ may be absorbed in and simply lead to a irrelevant 
change in frequency that will not affect the coefficients Fk{y) / Fk{Xk) of the 
wave function (10) which are independant of which eigenvalue is normalised. 
Also the zero in (12a) contributes nothing to the wave function and the theorem 
counts degenerate eigenvalues only once. 

The 36 X 36 matrix ^ decomposes into two bisymmetric diagonal 18 x 18 matrices 
of which only one gives the required eigenvalues (12a). By symmetry we choose 
//gg = ji as the matrix element corresponding to the ground state and then 
using (10) find a unitary wave function 

256V'9g = exp(i|6'){74 sin 6^ + 38 sin f 6* + 50 sin |6l + 25 sin 1461}, (13) 

dil)gg / de\g^^ = ji 

Again there is a high frequency spin orbital sin \A9 associated with the 5 orbiting 
electrons with 25 exchange correlations yielding an exchange correlation energy 
per electron of —24.6539/25 — —0.9862 hartrees using the results of |^. 

Finally the coherent electron coefficient corresponding to a zero spin state of C 
is 

720C6 = 7*^ + 1407"* + 4144A2 + 14400 (14) 

and ^ = 1 — 6*6/84, has the degenerate eigenvalues 

{12;|;l;if;0} 

belonging to a 24 x 24 bisymmetric matrix with ^12,12 = 211/84 labelling the 
ground state in which (10) yields a wave function 

256V'i2,i2 = sin 6I{39 sin if 6* + 128 sin 6* + 26.5 sin |6I + 37.5 sin 126*} (15) 



In this example which is not a metal, it would seem that there is an elec- 
tron exchange shared by the adjacent state sin |0 such that the total number 
of exhange correlations associated with the outer orbitals remains constant at 
26.5 + 37.5 = 64. If we assume that the largest negative exchange correlation 
energy is favoured then the exchange correlation energy per electron will be 
-37.8450/38 = -0.9959 hartrees. 

We will conclude this section with a comparison of the exact ground state cor- 
relation exchange energies with those derived by Gell-Mann and Brueckner for 
metals ||j who give the formula 



0.0622 In r, +0.094 hartrees 



_ _1 /" o.916 
'ff ~ 2 \ r^ 

Where Tj is the inter-electron spacing in units of Bohr radius. Thus if we take 
volume occupied by one electron to be r^ then the inverse density is A''"^ = r^ 
in these units, where N = A is the number of electrons, and (16) becomes 

eg = -i (0.916iVi/3 + 0.02073 In ]V + 0.084) hartrees 

Here the constant term has been slightly modified so that Cg = — 2 hartrees 
when N = 1 which is the known value for hydrogen. 

Table II compares eg values derived from (16a) with those calculated in sections 
2,3. 

Table II 
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Atom 


Exact eg 


eg by (16a) 


H 


-0.5 


-0.5 


Be 


-0.9778 


-0.9746 


B 


-0.9862 


-1.0476 


C 


-0.9959 


-1.1113 



4 Conclusion 



According to Biedenharn and Louck |g|. Section 7.7 this work is, in fact, a 
density functional formaHsm because the wave function X^^a^a is generated 
by a left ideal and (10) has the form of a general density matrix because if 
AA;Ffc(/i)//iFfe(Afe) is idempotent then so is Fk{^J,)/Fk{Xk) since Fk{n) vanishes 
unless ^ — k. The use of density matrices to describe canonical ensembles 
has been detailed by Prigogine [Q. In fact the idempotent matrices are Slater 
determinants and principal minors which describe pure or coherent states, so 
one is looking at nett many particle behaviour instead of attempting to build 
this up from the interaction of individual particles. Specifically we have the 
partition 



{A-X) + X = A 

into the two states of a canonical ensemble namely; {A — A) electrons with pos- 
itive spin and A electrons with negative spin as it is possible that very Hght 
elements such as Li may only be approximated. The intimate connection be- 
tween statistics and idempotency was actually recognised long ago by Eddington 
13 , Section 66. 
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